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The behaviour of magnetic perturbations of an initially potential three-dimensional equilibrium 
magnetic null point are investigated. The basic components which constitute a typical disturbance 
are taken to be rotations and shears, in line with previous work. The spine and fan of the null 
point (the field lines which asymptotically approach or recede from the null) are subjected to such 
rotational and shear perturbations, using three-dimensional magnetohydrodynamic simulations. It 
is found that rotations of the fan plane and about the spine lead to current sheets which are 
spatially diffuse in at least one direction, and form in the locations of the spine and fan. However, 
shearing perturbations lead to 3D-localised current sheets focused at the null point itself. In addition, 
rotations are associated with a growth of current parallel to the spine, driving rotational flows and 
a type of rotational reconnection. Shears, on the other hand, cause a current through the null which 
is parallel to the fan plane, and are associated with stagnation-type flows and field line reconnection 
across both the spine and fan. The importance of the parallel electric field, and its meaning as a 
reconnection rate, are discussed. 
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I. INTRODUCTION 

Many important physical phenomena in astrophysical 
plasmas are powered by magnetic reconnection. The 
locations at which reconnection is likely to occur in a 
complex three-dimensional magnetic field are those re- 
gions where strong currents (possibly singular in the ideal 
regime) may develop. However, determining what these 
locations might be is a non-trivial problem. 

There are a growing number of theories and pieces of 
evidence to suggest that 3D null points and separators 
(magnetic field lines which join two such nulls) may be 
sites of preferential current growth, in both the Solar 
corona and the Earth's magnetosphere. The field topol- 
ogy in the vicinity of such a 3D null point is defined by 
the field lines which asymptotically approach (or recede 
from) the null. These fall into two categories. A single 
pair of field lines approach (recede from) the null from op- 
posite directions, and are termed the 'spine' (or 7-line), 
while an infinite family of field lines recede from (ap- 
proach) the null in a surface called the fan (or £-) plane. 
The spine and fan of a given null may be determined by 
examining the linear field topology near the null, defined 
by the equation 

B = Mr, 

where the matrix M. is given by the Jacobian of B [e.g. 
Hoj . The eigenvectors of M. (whose corresponding eigen- 
values sum to zero since V • B — 0) now define the spine 
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and fan. The two eigenvectors whose eigenvalues are of 
like sign (or whose real parts have like sign) lie in the fan 
plane, while the third points along the spine. The fan 
surface is a separatrix surface of the magnetic field, sep- 
arating unique topological regions. While the spine does 
not separate topological regions (being only a line), it is 
nonetheless an important geometrical feature of the field. 
This is because, firstly, magnetic field lines converge on 
(or diverge from) it, and secondly, it is a field line along 
which disturbances may be channelled towards the null 
(as we shall see later). 

3D null points are predicted to be present in abundance 
in the solar corona. It should however be noted that 
many recent models of the magnetic field above the so- 
lar surface are based upon the Magnetic Charge Topology 
approach, extrapolating a potential field from point mag- 
netic sources in the photosphere. This is clearly a great 
simplification, and many (or all) of the 'photospheric null 
points' in such models will not be present when these 
idealised flux sources are replaced by more realistic finite 
flux patches. However, various approaches also predict 
the presence of 3D null points up in the corona, with an 
average of between approximately 7 and 15 coronal nulls 
ex pec ted for every 100 photospheric flux concentrations 
P.ll8L[3l|. The presence of these coronal nulls is expected 
to be more robust to the method of field extrapolation 
0, Q • The separatrix surfaces and separators associated 
with these coronal null points are thought to be likely 
sites of coronal heating and reconnection [l?], H3] ■ F ur_ 
thcrmorc, there is observational evidence that reconnec- 
tion at a 3D null point (both fan-type and spine-type 
reconnection) may be at work in some solar flares [9(. 
Closer to home, separator reconnection is thought to oc- 
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cur on the dayside of the Earth's magnetosphere [e.g.l33|. 
In addition, there has been a recent in situ observation 
by the Cluster spacecraft [33[ of a 3D magnetic null, at 
which it is proposed that reconnection is occurring, lo- 
cated within the current sheet in the Earth's magnetotail. 

Although it is now realised that 3D nulls are of great 
importance for reconnection in realistic 3D geometries, 
what is still lacking is a clear picture of what reconnec- 
tion processes at such nulls look like. In particular, a 
description of the sorts of physical signatures expected 
(current sheets, plasma flows etc.), in order to lead the 
analysis of new observations, is of great importance. Our 
aim in this paper is to go some way to providing such a 
picture. 

While the relationship between reconnection at a sep- 
arator and a single 3D null point is not well known, it is 
clear that the two should be linked in some way. What 
is known is that both areprone to collapse in response 
to external motions [e.g. [H, [22[. That is, a 3D null 
point is a 'weak point' of the magnetic field to which 
disturbances are attracted [UEl- Furthermore, kine- 
matic considerations suggest that null points and sepa- 
rators are both locations where singularities may form in 
ideal MHD [e.g. [28[. The attraction of disturbances 
to 2D nulls is also well documented [e.g. M, [T^|. 

The local behaviour of a perturbed 3D magnetic null 
point has been examined in the linear regime for the cold, 
resistive MHD equations by [29[ . They employed a modal 
decomposition, and determined that only m = and 
m = 1 modes can lead to currents at the null point, 
where m is the azimuthal wavenumber. In addition, they 
found that while m = modes are attracted to the spine 
and fan of the null, m = 1 type perturbations tend to 
focus in towards the null itself. This was demonstrated 
in the linear regime with 2D simulations in the rz-plane. 
One of the major aims of the present work is to ascertain 
whether such behaviour for the evolution of disturbances 
is found in the full MHD regime, using 3D simulations. 

In addition to the above, we will investigate which cur- 
rent components develop at the null in response to dif- 
ferent perturbations, and the implications which this has 
for the plasma flow and field line behaviour. In each 
case we consider, the plasma is initially at rest, and so 
with no flow through our boundaries, we do not com- 
pare with strongly driven flux pile- up models [e.g. It 
is anticipated that different current orientations at the 
null will lead to very different behaviour 0, [25| . One of 
the present simulations is very similar to that described 
in the work of [l2l ]. and so is only briefly summarised in 
what follows. The other simulations combine with this 
to create a complete picture of the subject at hand. 

In Section |TT] we describe the numerical scheme and the 
initial conditions used. In Sections |TTT] [VI] we describe 
the results of our simulations for various different types 
of perturbations, and in Section [VIII we give a summary 
and conclusions. 



II. NUMERICAL SCHEME AND SIMULATION 
SETUP 

A. Numerical Scheme 

The numerical scheme employed in the simu- 
lations which follow is described briefly below 
(further details may be found in [21[ and at 
|http://www. astro. ku.dk/~kg| ). We solve the three- 
dimensional resistive MHD equations in the form 
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where B is the magnetic field, E the electric field, v the 
plasma velocity, 77 the resistivity, J the electric current, 
p the density, r the viscous stress tensor, P the pressure, 
e the internal energy, Q V i SC the viscous dissipation and 
Qj the Joule dissipation. An ideal gas is assumed, and 
hence P = (7— 1) e = |e. 

The equations above have been non-dimensionalised by 
setting the magnetic permeability p,o = 1, and the gas 
constant (TV) equal to the mean molecular weight (M). 
The result is that, for a cubic domain of unit size, if \p\ = 
|B| = 1, then time is measured in units of the Alfven 
travel time across the domain (T4 = L^/pp a /Bo, where 
L is the size of the domain, and po and Bq are typical 
values of the density and magnetic field respectively). 

The equations are solved on staggered meshes; in this 
way the required MHD conservation laws are automat- 
ically satisfied. Spatial derivatives are evaluated using 
a sixth-order-accurate finite difference method. It is of- 
ten the case that, due to the staggered mesh, this value 
of the derivative is returned in exactly the position it is 
needed. When this is not the case, the values are cal- 
culated using a fifth-order interpolation method at the 
relevant position. A non-uniform mesh is used in the 
simulations, with higher resolution in the centre of the 
domain in each coordinate direction. In this way it is 
possible to better resolve behaviour in the vicinity of the 
null, and indeed across the spine and fan, while keep- 
ing the boundaries 'far away' to reduce their effect. A 
third-order predictor-corrector method is employed for 
time-stepping. All simulations are carried out on grids of 
128 3 or 256 3 resolution. 

Viscosity is handled using a combined second-order 
(constant v) and fourth-order method, which is capable 
of providing sufficient localised dissipation where neces- 
sary to handle the development of numerical instabili- 
ties. The result is that it is possible to achieve much 



lower effective values of v for a given numerical resolu- 
tion than with a constant-^ approach. This is achieved 
by having an enhanced viscosity at length scales close 
to the numerical resolution limit, effectively dissipating 
short wavelength disturbances, while leaving larger-scale 
structures nearly undamped. Viscosity defined in such a 
way is often termed 'hyper-viscosity'. Such an approach 
is used in 3D simulations to maximise the fraction of the 
3D domain that has an 'ideal' behaviour. Two different 
models are used for the resistivity; either a traditional 
constant-^ (second-order) model or a hyper-resistivity 
model (labelled r\h in the text) similar in approach as 
the viscosity model. Comparisons between identical sim- 
ulation runs using the two resistivity approaches show 
the same general evolution, with the main difference be- 
ing that the hyper-resistivity runs show more spacially 
localised structures. This is demonstrated explicity for 
one particular case below (see Section W} . 

In a number of the simulation runs, 'trace particles' are 
used in order to track the motion of field lines in time. 
These points are chosen at t — 0, and are then followed in 
time throughout the simulation by integrating the plasma 
flow field. These 'fluid elements' are then used to define 
field line footpoints. 

Finally, the boundary conditions are closed in all three 
directions. In addition, as we aim to study only the ini- 
tial localisation of perturbations at the null, a boundary 
damping zone is included, in order to limit the reflec- 
tion of waves back into the domain from the boundary. 
Within this region, a fraction of the kinetic energy is re- 
moved per unit time, which, when chosen appropriately 
for the wave speed, may effectively damp a large por- 
tion of the wave energy. The effect of the boundaries 
on the dynamics of the null point is, on the other hand, 
negligible, as the simulations are all terminated before 
the perturbations have time to reflect off them and reach 
back to the region around the null again. 



B. Initial setup 

The simulation is set up as follows. The initial 'back- 
ground' magnetic field is potential, and describes a 3D 
magnetic null point located at the origin, specifically 
B = Bo (x, y, —2z). The plasma is initially at rest, with 
p = 1, e = 5/3/2 everywhere, where (3 determines the 
plasma-/?. The system is then disturbed by perturb- 
ing the magnetic field. We choose to add a magnetic 
field instead of, say, imposing some velocity on a line- 
tied boundary, so as to examine the evolution of the null 
point, and the reconnective behaviour, in an undriven sit- 
uation. We focus on the behaviour of the transient pulse 
as it moves towards the null. Note that the fact that 
there is no flow through our boundaries means that we 
do not compare with strongly driven flux pile-up models. 

In order for the null point to be affected, the pertur- 
bation must disturb either the spine or fan of the null, 
otherwise the disturbance will simply propagate back and 
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FIG. 1: Schematic view of isosurfaces of the different pertur- 
bation magnetic fields, with arrows indicating the field orien- 
tation, (a) rotation within the fan plane, (b) shear in the fan 
plane (c) rotation of the spine, in the same or opposite sense 
above and below the fan, and (d) shear of the spine. In each 
case the vertical line is the (undisturbed) spine of the null, 
and the square is the fan. 

forth along the associated magnetic field lines, bouncing 
between the boundaries. The perturbation magnetic field 
must of course be divergence- free, and in each of the cases 
discussed takes the general form 

6, = b R x cos(fli) exp (~ {Rl -J lo)2 - , 

where bi and bj are the two components of the perturba- 
tion magnetic field, £ is the third spatial coordinate, and 
bo, -Rio, Co j a>h a nd b^ are constants. R\ and Q\ are de- 
fined depending on the orientation of the perturbation. 
The perturbation is localised by the exponential terms 
within either a linear tube or a torus. 

Five different types of perturbation will be considered. 
These are described briefly below, and are illustrated in 
Figures Q] and [H Figure [1] shows isosurfaces of the dis- 
turbance field magnitude, with the arrows showing the 
direction of the field. Figure plots a selection of field 
lines from around the null, showing the effect of the dis- 
turbances on the the magnetic field. The first type of 
perturbation corresponds to the m = disturbances of 
[29j | . that is, rotational motions. We consider two cases, 
in the first of which (Section IIII[) the rotation is con- 
centrated in the fan plane (but away from the null, see 
Figures [Ha), Ufa)). One may also perform a rotation 
about the spine, but away from the fan plane, either in 
the same or opposite sense above and below the fan (see 
Figures [ljc), (He)). This has been investigated in detail 
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FIG. 2: Sample magnetic field lines at t — for the different 
perturbations, (a)-(d) as in Figure Q] 

by [Hj] i an d is discussed briefly in Section ITVl Alterna- 
tively, the null point separatrices may be perturbed by 
applying some shear (corresponding to m = 1), to ei- 
ther the fan (Section [V] Figures [T^b), HJb)) or the spine 
(Section IVI1 Figures [ljd), HJd)). Any more generic per- 
turbation may be made up from a combination of such 
rotations and shears. 



III. ROTATION IN THE FAN PLANE 

For a rotation in the fan plane, the perturbation mag- 
netic field lies in the rry-planc, such that {bi,bj} — 
{bx,b y } and £ = z in Equation l[7|). and in addition we 
take i?i = \J x 1 + y 2 , 6i = tan _1 (?//x), and Co = 0. 
Thus the current is initially concentrated in a torus whose 
toroidal axis lies in the fan plane at a radius Rio from 
the origin (null), see Figures Ufa) and a). The domain 
size is chosen to be 1.5 x 1.5 x 2.5, in order to limit the 
effect of the upper/lower boundaries in z. We use the 
hyper-resistive model, r\ = r\h, and the characteristic pa- 
rameters for the experiment are; Bq = 1, bo = 0.1, Rio = 
0.16, (3 — 0.01, ah = bh = 0.06, giving a travel time to the 
top (z) boundary of order 1.8 and to the xy-boundaries 
of 1.5 in code units, for the main body of the disturbance. 



A. Current density evolution 

Figure [3] shows the time evolution of the current mod- 
ulus in a plane through the null and spine axis (cho- 
sen to be the y = plane, although we have cylindrical 
symmetry away from the boundaries) . It is evident from 
the image that pulses propagating both inwards and out- 
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FIG. 3: For a rotation of the fan; (a) current modulus in 
the y = plane for [x,z] = [±0.75, ±0.75], at times marked 
by asterisks in (b) (t = 0,0.72, 1.44,2.88). (b) Evolution of 
the maximum value of each current component (J x and J y 
dashed, J z solid, black) and of the integrated parallel electric 
field along the spine (grey). Domain size is 1.5 x 1.5 x 2.5 
and parameters used are Bo = 1, 6o = 0.1, -Rio = 0.16, j3 = 
0.01, a;, = 0.06, bh ~ 0.06 and r\ = r\h- The dot-dashed line 
shows J z evolution in a run with the same parameters, but at 
256 3 resolution. 



wards develop, as one would expect. In the following, we 
refer to the disturbances which propagate towards and 
away from the null point as the 'ingoing' and 'outgoing' 
pulses, respectively. Note that here, as in the simulations 
which follow, both the ingoing pulse and outgoing pulse 
have a somewhat complicated structure. Specifically, it 
appears from Figure G^a) that there are two wavefronts 
which localise towards the null, although in fact these 
two wavefronts (in |J|) correspond to current concentra- 
tions of opposite sign (in J z ), which demark strong field 
gradients at the 'back end' and 'front end' of the single 
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ingoing 'pulse'. The magnetic field gradients must neces- 
sarily pass through zero at the centre of this pulse, hence 
the appearance of two strong bands in |J|. The ingo- 
ing pulse propagates along the 'background' field lines, 
and concentrates in a current tube centred on the spine. 
There is, however, no preferential attraction of the cur- 
rent to the null itself. It is interesting that even though 
the field structure is hyperbolic, the current localises very 
uniformly in z, i.e. the wavefront is very close to vertical, 
as would be expected for an Alfven wave. Note finally 
that in the last image there are some unavoidable reflec- 
tions of the disturbance from the upper boundary back 
along the background field lines. We have checked that 
the boundaries play no significant role in the qualitative 
or quatitative evolution by repeating the simulation in a 
domain of dimensions 1.5 x 1.5 x 4.5 (which increases the 
perturbation travel time to and from the z-boundaries to 
greater than 4 in code units). 

Further insight may be gained by examining the evo- 
lution of the different components of the current. The 
maximum value of each component is plotted in Fig- 
ure Eljb) (we take the maximum within the region 
[±0.25, ±0.25, ±0.42] in each direction to exclude bound- 
ary effects) . The figure shows that as the current localises 
towards the spine, it is J 2 , i.e. the current parallel to the 
spine, which is significantly amplified, while J x and J y 
are not. This is to be expected, as it is dx and dy which 
are decreasing as the disturbance is squeezed in towards 
the spine. The current eventually reaches a maximum 
as it localises. This is due to a combination of the finite 
numerical resolution and the imposed resistivity model. 
For constant 77, increasing the numerical resolution (TV) 
will eventually result in a behaviour of the solution that 
is independent of N. On the other hand, the current in 
the hyper-resistive case will continue to increase with N, 
with the physical structure having a length scale com- 
parable to the numerical resolution. This can be seen 
by examining the growth of the dominant (z) component 
of J for a run with double the resolution, shown by the 
dot-dashed line in Figure \3j[b). The early evolution is 
very similar, though a higher current peak is eventually 
achieved after a longer period of localisation. 

In a truly ideal evolution, it is expected that the cur- 
rent would increase indefinitely in time, although more 
and more of the energy associated with the disturbance 
would 'escape' down the spine [c.f. [2(|. However, in a 
true physical case, no matter how small the resistivity, it 
will always become important eventually, once the cur- 
rent becomes sufficiently intense, and as a result some 
energy can be dissipated. This will occur all along the 
spine. 

One very interesting question is how the maximum cur- 
rent depends on the resistivity. It is impossible here for 
us to use realistic resistive parameters appropriate for, 
e.g. the Solar corona, where the resistivity is typically of 
the order 10 -14 . One crucial feature of any reconnection 
model is therefore whether the reconnection rate scales 
as some negative power of 77 — if so then it will typically 



correspond to 'slow' reconnective behaviour in a realistic 
plasma. It is our intention to investigate this dependence 
in a subsequent paper in this series. It has furthermore 
recently been demonstrated @, [HJ that the plasma pres- 
sure (i.e. the value of the plasma-/?) can have a profound 
effect on current scalings at 2D and 3D null points. This 
will also be investigated further in a later paper. 

At a first glance, the type of tubular current structure 
described above is reminiscent of those found in the in- 
compressible 'spine reconnection' solutions first described 
by [5j (in their case there exist two tubes in close proxim- 
ity due to the assumed symmetry). In a general incom- 
pressible time-dependent case, a single tube of current is 
stretched out along the spine of the null, with the current 
being purely azimuthal within the tube (that is, directed 
parallel to the fan plane), rather than axial [23]. This is 
clearly rather different from the situation we have here, 
where the current is very close to being parallel to the 
spine. In addition, the incompressible spine reconnec- 
tion solutions are associated with reconnection of field 
lines across the spine and fan (advection across the fan 
and diffusion across the spine), which we shall see below 
does not occur here. 

One final important characteristic of the current evo- 
lution observed in our simulation is a non-linear coupling 
which may occur when the perturbation magnetic field 
is strong compared with the background field. As ob- 
served by the main Alfvenic disturbance may also 
couple to a fast-mode wave which is attracted to the null 
point itself. This is, however, a fairly weak effect. [l2j 
found that the coupling only existed in cases where the 
boundary driving was ramped sharply up from zero, and 
similarly we find that the strength of (and gradients in) 
the disturbance field must be very strong in order that 
the effect can be seen at all. 



B. Plasma flow and field line behaviour 

An examination of the plasma flow induced by the 
perturbation is of interest in itself, as well as in help- 
ing to determine what type of reconnective behaviour 
might result once a current sheet forms. It is perhaps 
not surprising that the rotational perturbation induces 
rotational plasma flows centred on the spine, in the xy- 
plane, as shown in Figure Ufa). The rotation is present 
within the 'envelope' of magnetic flux which was magnet- 
ically connected to the initial perturbation. Importantly, 
there is no plasma flow across either the spine or fan of 
the null. In fact, the direction of the plasma rotation in 
Figure HJa) is somewhat complicated by reflections from 
the z-boundaries. Initially, the ingoing pulse generates 
clockwise flow (viewed from the z-boundaries) and the 
outgoing pulse anti-clockwise. In addition, once the in- 
going pulse is reflected off the z-boundaries, it generates 
an additional region of anti-clockwise flow, which prop- 
agates in the negative z-direction. This occurs first at 
large radius within the flux envelope since this is where 



6 



-0.5 




FIG. 4: Plasma flow in the plane z ~ 0.05, showing rotation 
centred on the spine (centre), for perturbations (a) rotation 
in the fan and (b) rotation about the spine. Pattern in other 
planes z = zo is similar, with only the r-localisation varying 
due to the hyperbolic nature of the flux enevelope affected by 
the perturbation. The background shading shows the current 
modulus in the same plane in each case. Each image is at 
the time of maximum current, and the parameters are as in 
FigureEJfor (a), and for (b) B = 1, b = ±0.1, Co = ±0.2, /3 = 
0.01, a h =0.06,fe h =0.06, 77 = rj h . 



the reflection first occurs (see Figure[3]). Thus, present in 
Figure \M&) are the flow resulting from the ingoing pulse, 
and its reflection from the z-boundary. 

As a result of this plasma flow, it seems unlikely 
that any reconnection involves field lines being advected 
across the separatrices, but rather a rotational mismatch- 
ing would be expected. This is predicted by kinematic 
studies [25[ when the current is parallel to the spine of 



the null point. As shown in Figure [5l this is indeed the 
case. The plotted field lines are traced from two sets of 
trace-particles, which are initially connected (see Figure 
Uta)). The first set of particles (from which the grey (blue 
online) field lines are traced) remain approximately in 
the ideal region, i.e. they are located far away from the 
null (near the fan plane) where currents remain weak. 
These field lines therefore show approximately the ideal 
behaviour, which is followed everywhere except close to 
the spine. (When the field lines pass close to the spine 
in the figure it is not possible to see the individual lines, 
due to the converging field structure.) The other set of 
particles is located close to the spine, such that they are 
eventually engulfed by the localising current. It can be 
seen that the corresponding field lines (black, red online) 
continually change their connections [see [26| in a rota- 
tional fashion as the current localises. Note that in the 
simulation pictured, a higher value of bo, and a constant 
resistivity (77 = 0.002, implying a Lundquist number of 
order 500) have been used for illustrative purposes, since 
the diffusive region is larger for constant 77 than with the 
hyper-resistivity, and the effect of the diffusion is there- 
fore easier to visualise. 



C. Parallel electric field 

It is well-known that a crucial indicator of 3D magnetic 
reconnection is the presence of an electric field compo- 
nent parallel to the magnetic field, Ei\ [e.g. 30]. Examin- 
ing an isosurface of E\\ , it is clear that non-ideal effects 
become important basically uniformly all along the spine 
(see Figure |6j taken a little before current maximum to 
limit appearance of boundary effects). Thus the recon- 
nection seems to be associated with a spatially diffuse (at 
least in one dimension) region. 



24j showed that in the case of an isolated 3D null with 



current parallel to the spine, the rate of rotational flux 
mis-matching can be quantified by calculating the inte- 
grated parallel electric field (from one end of the diffusion 
region to the other) along the spine of the null (see also 
[l5j). Care must be taken, however, in comparing our 
simulations with this kinematic result. The principal rea- 
son for this is that the result relies upon the assumption 
that the diffusion region is bounded in z, encompassing 
the null. Nonetheless, we find that $ s = J En (along 
the whole length of the spine within the domain) does 
indeed show a clear peak in time, indicating a maximum 
in the reconnection rate at the null, which occurs once 
the perturbation reaches the spine. That is, when the 
field gradients built up sufficiently that the current sheet 
begins to diffuse onto the spine, a parallel electric field 
can be seen to develop along it (see Figure [3]). While the 
numbers we obtain here (Figure [3fb) ) for En are depen- 
dent on the resistivity model and/or value, the qualita- 
tive behaviour is not. 

It is interesting to observe the close temporal corre- 
lation between the peak reconnection rate (3> s ) at the 



FIG. 5: Rotational slippage between two sets of magnetic field lines. The black loops enclose the two sets of fluid elements 
from which the field lines are traced. For parameters Bo = 1, 60 = 3, Rio = 0.18, j3 = O.Of , an = 0.06, bh = 0.06 and constant 
r; = 0.002. Again the total domain size is 1.5 x 1.5 x 2.5. Times for the plots are t = 0,0.2, 0.9. 
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FIG. 6: Isosurface of E« resulting from a rotation in the fan 
plane, at 50% of its maximum, at t = 2.16 and for the same 
parameters as Figure |3] 



null and the peak current. As mentioned previously, in 
the ideal limit we would expect J to grow indefinitely, 
due to the structure of the magnetic field in the vicin- 
ity of the null. Thus a current peak is also an indication 
(for our transient perturbations) that non-ideal processes 
have become important, allowing the stress in the field 
(here twist) to begin to significantly dissipate. With this 
type of perturbation, Eu will also increase as the pulse 
approaches the spine due to the geometry of the mag- 



netic field, which becomes increasingly .E^ -dominated as 
one approaches the spine line (at a given height, z). 

IV. ROTATION ABOUT THE SPINE 

For a rotation about the spine we may take a per- 
turbation magnetic field of similar form to that used in 
the previous section. However, here we take £0 non-zero 
and Rio to be zero. For perturbation above and be- 
low the fan plane, we superpose two disturbances of this 
form, with £q of opposite signs (Figures [TJc) and[^c)). 
Whether the rotation is in the same or opposite sense on 
either side of the fan is determined by the relative signs 
of bo. We take B = 1, b = ±0.1, = 0.01, Co = ±0.2, 
dh = bh = 0.06. A thorough examination of the devel- 
opment of these types of perturbations has been carried 
out in a similar simulation, described by [l2l ]. The one 
major difference between the two simulations is that they 
applied a driving velocity at the (line-tied) boundaries, 
whereas we perturb the magnetic field within the domain. 
We therefore only summarise the results here. 

This time, as the perturbation evolves, it generates 
a current which spreads out in the fan plane as it ap- 
proaches the null (see Figure [7ja)). However, the most 
important aspect of the current development is that there 
is once again no tendency for a focusing of the localisa- 
tion towards the null itself. Regardless of whether the 
rotation has the same or opposite sense, there are inward 
and outward propagating disturbances, and in the ingo- 
ing disturbance 3 xy is amplified (Figure [3b, d)). The 
outgoing pulse, which localises towards the spine in the 
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same way as was seen in Section IIII1 shows a growth of 
J z (Figure HJc)). However, this time that current never 
reaches the null since the fan is not perturbed. There is 
also a relatively weak pulse of J z which does propagate 
towards the fan plane and null (JZt c )), as seen by [29j . 
however, this is dominated by the strong J xy further out 
in the fan. 

The fact that different current components are mag- 
nified in the ingoing and outgoing pulses is likely just a 
result of the geometry of the background field. In par- 
ticular, moving inwards from a given point on the spine, 
d/dx,d/dy decrease, while d/dz increases sharply as the 
fan is approached. Thus the current maximum occurs as 
the ingoing pulse steepens towards the fan. If the resis- 
tivity were zero, a fan current sheet (with infinite cur- 
rent density but of zero thickness) would be asymptoti- 
cally approached. However, in our simulation, sufficiently 
steep gradients develop before the disturbances reach the 
fan, such that diffusion becomes important. This is one 
difference between our simulation and the driven case of 
[l2| , where the continual twisting at the boundaries drove 
the current right into the fan plane. 

The plasma flow in this case is, as expected, of a rota- 
tional nature, with similar form to that found for rotation 
in the fan plane (and sense of rotation governed by the 
initial condition). A plot of the plasma velocity in the 
xy— plane demonstrating this is shown in Figure E^b). 
Furthermore, the field line behaviour is found to be very 
similar to that shown in Figure [5l 

As the disturbances here travel basically along B and 
there is no flow and no magnetic connection between the 
two sides of the fan plane (z > 0, z < 0), there is essen- 
tially no difference for the case where the driving has the 
same sign on each side of the fan, except the sense of ro- 
tation. In addition, since the current is not driven right 
into the fan as in [I2[ , there is no issue of the currents 
reinforcing or cancelling there as they found. 



V. SHEAR OF THE FAN PLANE 

The behaviour of shearing-type perturbations is very 
different to rotational ones. We first consider the case 
of a shear of the fan plane. The perturbation again 
takes the form described by Equation ([7]), where this 
time {bi,bj} = {b x ,b z } and £ = y. In addition we 

take R\ = \/(x — xq) 2 + z 2 , 9\ = tawT 1 ( — — ) and 

Rio = Co = 0. This corresponds to a linear tube (of 
finite length) of azimuthal magnetic flux, whose axis lies 
along x = xq, z = 0. We superimpose two such pertur- 
bations on our background null, with opposite signs of 
xq, as shown in Figure [ljb). The effect on the fan field 
lines is as if they have been 'plucked', in one direction 
in some region of the fan, and in the opposite direction 
in an opposite region (Figure EJb)). We take (except 
where stated) B = 1, b = 0.1, (3 = 0.01, x = ±0.16, 
a h = 0.06, b h = 0.2. 
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FIG. 7: For a rotation about the spine of opposite sense above 
and below the fan; (a) shaded images showing |J| in the y = 
plane, for [x,z] = [±0.45, ±0.45], at times shown by the 
asterisks in (d) (t = 0,0.45,0.9,1.35). (b) J x and (c) J z at 
t = 1.35 (time of maximum current), (d) Evolution of the 
maximum values of each current component in the inner half 
(in each direction) of the domain: J x dotted, J y dashed, J z 
solid line. For domain size 1.5 x 1.5 x 1.5 and parameters as 
in Figure H^b) . 



A. Current and plasma flow 

The crucial characteristic of the current development 
for a shear of the fan is that, while the disturbance 
propagates to some extent along the background field 
lines, there is additionally a strong focusing of the cur- 
rent (across the field) towards the null point itself (see 
Figure HJa)). Note that again there is a fair degree of 
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FIG. 8: For a shear of the fan; (a) shaded images showing 
J in the y = plane, for [x,z] = [±0.45, ±0.45], at times 
marked by asterisks in (b) it = 0,0.38, 1.32, 2.45). (b) Evo- 
lution of the maximum values of each current component in 
the inner half (in each direction) of the domain: J x dotted, 
J y dashed, J z solid black line. The grey line plots the in- 
tegrated parallel electric field along the y-axis. For domain 
size 1.5 x 1.5 x 1.5 and parameters Bo = 1, bo = 0.1, xo = 
±0.16, (3 = 0.01, o h = 0.06, b h = 0.2, r) = r\ h . The dot-dashed 
line shows J y evolution in a run with the same parameters, 
but with 77 = 5 x 10~ 5 , constant 



structure to the localising current. This is because the 
disturbance was initiated by a divergence-free magnetic 
disturbance, in the form of a tube of magnetic field (see 
Figure |Tfb)) . Thus in fact three current concentrations, 
demarking field gradients at the front, middle and back 
end of the pulse are present. The components at the front 
and back end are necessarily of opposite sign (in J y ) to 
the strongest concentration in the middle, and the gra- 
dients of course go through zero between these regions, 
hence the appearance of three 'pulses' in |J|. While the 
two outer regions initially join up to encircle the inner 



concentration, this appearance is eventually lost due to 
the stretching of the structure in z. 

It is evident from Figure [Htb) that the disturbance 
which focuses towards the null point is dominated by 
J y , that is, the current orthogonal to the plane of the 
shear. Once again it is the resistivity which limits the 
growth of the current. In an ideal situation we would 
expect the current to continue to grow indefinitely, fo- 
cusing at the null, where a good portion of the energy of 
the disturbance would be deposited, once resistive effects 
eventually become important. 

As an explicit demonstration of the effect of our hyper- 
resistivity model, we re-run the simulation with r\ con- 
stant, taking 77 = 5 x 10~ 5 (in order to resolve all struc- 
tures sufficiently). We find that the diffusion causes the 
peak current to be much lower, and occur much ear- 
lier when the perturbation is far less localised (see the 
dot-dashed line in Figure [8jb)). Qualitatively the cur- 
rent structure behaves very similarly — the same images 
as in Figure 0(a) look almost identical (when scaled by 
their individual maxima), except that in the final image 
(now long after current maximum), the current is not so 
strongly peaked near z = as much of the perturbation 
has diffused away. 

The presence of strong current at the null point di- 
rected parallel to the fan is expected to indicate the 
presence of more traditional reconnection-type flows. In 
the simplified kinematic model of 25], a stagnation-point 
flow, centred on the null, is present in the plane per- 
pendicular to the shear. This flow transports magnetic 
flux across the spine and fan. Similarly here, at the 
time of maximum current (see Figure [5]( a)), we indeed 
find strong plasma flow across the fan plane (originally 
z = 0), concentrated around the areas of maximum cur- 
rent. Away from the current concentrations as well, there 
is flow across the original locations of the spine and fan. 
However, it is not straightforward to determine whether 
magnetic reconnection is occurring at the null. This is be- 
cause the flow may act to transport magnetic flux across 
the fan plane (magnetic reconnection), or it may simply 
act to advect the spine and fan in an ideal sense. It is 
therefore again of interest to examine the development of 

E \\- 



B. Eu and reconnective behaviour 

The presence of a parallel electric field denotes a break- 
down of ideal behaviour, and so is crucial in determining 
whether magnetic reconnection occurs. We expect E\\ 
to be greatest in regions of high current, where the field 
lines lie in the direction of the current (since E\\ = r]J\\). 
Examining an isosurface of Eu at current maximum (Fig- 
ure [TU)) . we see that in fact the highest concentrations are 
located (as expected) at regions of strong J, a little way 
from the y-axis (due to the structure of the initial distur- 
bance) . The most striking thing we see from Figure [TU] 
is that in this case the regions in which Eu develops are 



10 




-0.16 



-0.28 



0.28 



0.6 



-0.6 



-0.6 



x 



0.6 



FIG. 9: (a) The plane y = 0, perpendicular to the shear, 
at time of maximum current (t = 2.45), for a shear of the 
fan. The arrows show the plasma flow, and the background 
shading the current modulus. Parameters as in Figure [8] (b) 
As (a), but for spine shear run (t = 1.23), for domain size 
1. 5x1. 5x1. 5 and parameters Bo = 1, bo = 0.1, zq = ±0.2, (3 = 
0.01, a h = 0.05, b h = 0.05, t] = rj h . 



highly spatially localised in all directions. Hence recon- 
nection processes are very local to the null itself. While 
this localisation is not quite so extreme in runs with con- 
stant resistivity, the basic structure is the same. 

For an isolated null with fan-aligned current, sur- 
rounded by a localised non-ideal region, the integral of 
JSii along the magnetic field line in the fan plane which is 
directed along J can be shown to give an exact measure of 
the rate of flux transfer across the fan, where the integral 




FIG. 10: For a shear of the fan plane, (a) isosurface of E\\ at 
25% of its maximum and (b) surface of E\\ in the z = plane 
(interpolated onto a uniform grid), both at t = 2.45 (current 
maximum), for the same parameters as in Figure [S] 



is taken from one side of the diffusion region (assumed 
spatially localised at the null) to the other [25|]. Thus, 
a non-zero value for <&/ = J Ends along a fan field line 
threading the diffusion region implies that magnetic flux 
is transferred across the separatrix (fan) plane, even in 
the ideal region, and thus demonstrates that reconnection 
is occurring at the null. Performing such an integration, 
$/ = J x _ z _ Q E\\ds, we find a parallel electric field does 
indeed develop, a strong indicator that some reconnective 
process is taking place. Note that here again the peaks 
of <f>/ and J max are closely correlated. 

In fact, when $y reaches its temporal maximum, the 
spatial maximum of E\\ is still some distance from the 
y-axis (see Figure [T0|) . In the fan plane, the maxima ac- 
tually occur on field lines which lie approximately along 
x = ±y. The reason for this is as follows. The electric 
field, like J, is directed largely in the y direction per- 
pendicular to the shear, and is uni-directional through 
the null. Thus E\\ w on the x-axis, and E\\ > for 
y > 0, Eu < for y < (see Figure HUfb)). If the per- 
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turbation were to extend azimuthally around the entire 
fan plane, we would expect the maximum and minimum 
of E\\ to be on the y-axis. Such an azimuthally sym- 
metric perturbation could be envisaged by considering 
a field localised within the torus of Figure [Ha) which 
is purely poloidal rather than toroidal, and which has a 
cos{6) multiplying factor (where 6 = tan -1 (y / x)) . How- 
ever, since the chosen shear has a finite azimuthal ex- 
tent, maxima and minima of E\\ are found off the y-axis. 
Therefore, by comparison with [25[, we expect a better 
measure of the rate of flux transfer across the fan to be 
given by $/ = j x=y ^ z=Zfan E\\ds, by symmetry, giving 
the maximal value of This quantity shows a similar 
qualitative evolution to that of plotted in Figure[5f b). 

The nature of the magnetic reconnection associated 
with En can again be examined by tracing field lines from 
trace particles which are initially magnetically connected. 
A typical evolution for an appropriately chosen set of field 
lines is shown in Figure QT] It is clear that the chosen 
plasma elements change their connections in such a way 
as to suggest a transport of magnetic flux across both 
the spine and fan of the null. Specifically, grey field lines 
(green online, traced from 'non-ideal' footpoints located 
near the fan plane, between the initial perturbation and 
the null) are transported across the fan plane (from above 
to below), whereas black field lines (red online, traced 
from 'ideal' plasma elements located far up the spine) 
are advected across the spine. 

VI. SHEARING THE SPINE 

To complete the picture we finally consider a pertur- 
bation which shears the spine of the null point. The 
disturbance takes the form given by Equation iff]), with 
{bi,bj} = {b x ,b z }, C = V, Ri = \Jx 2 + {z ~ Zq) 2 , 
9i = tan^ 1 ( £ ^ £a ) and R w = Co = 0. Again two such 
disturbance fields are superimposed, with opposite signs 
of z n (see Figure Did)). The effect on the spine field line 
is as if it has been locally 'plucked', in opposite direc- 
tions above and below the null (see Figure [UJd)). We 
take (except where stated) Bo = 1, bo = 0.1, f3 = 0.01, 
zo = ±0.2, a h = b h = 0.05. 

The resulting current evolution in the plane perpen- 
dicular to the shear is shown in Figure 1121 from which 
it can be seen that the inward-propagating component of 
the disturbance is again attracted towards the null point, 
as in the case of a shear of the fan. Although both J x and 
J y are magnified during this localisation, in fact J = J y 
at the null point itself. That is, the current which devel- 
ops at the null is again parallel to the fan, perpendicular 
to the plane of shear. 

As before, a plot of the plasma velocity in the plane 
perpendicular to the shear (see Figure E^b)), suggests 
flow across the spine and fan, although again it is hard 
to say whether the flow actually crosses the spine and 
fan, or merely ideally advects them. In addition, there is 
also some complicated (and relatively strong) flow across 



the spine within the region of the localising current con- 
centration. 

Some representative field lines are plotted in FigurefTBI 
Field lines traced from plasma elements initially located 
near the spine (black, green online) are transported back 
and forth across the spine, since these plasma elements 
are engulfed by the current concentration. The grey field 
lines (red online), however, are traced from plasma el- 
ements which stay forever in the 'ideal' region, far out 
along the fan plane, and are clearly advected across the 
separatrix surface. In fact, the behaviour of the null in re- 
sponse to the perturbation is similar to the case where the 
fan plane is sheared. In each case, the current is peaked 
at the null itself, and field lines reconnect in a manner 
similar to the traditional 2D picture, in the plane of the 
shear. 



VII. SUMMARY 

The behaviour of perturbations in the vicinity of a 
potential 3D magnetic null point has been investigated 
via resistive magnetohydrodynamic simulations. Distur- 
bances which affect the null point must deform the field 
in the vicinity of either the spine or the fan of the null. 
A typical such disturbance may be constructed using a 
combination of shears and rotations, and hence we con- 
sidered four basic possibilities — namely shear of the fan 
or spine, and rotation of the fan or about the spine. It 
should be pointed out that we consider a case with no 
initial 'background' flow, that is there is no flow through 
our boundaries, and therefore we do not compare with 
strongly driven flux pile-up models. 

We have found that rotational types of perturbations 
tend to lead to current accumulation in the vicinity of 
the field lines which asymptotically approach (or recede 
from) the null point, that is the spine and fan field lines. 
The disturbances behave in a way which is essentially 
Alfvenic, propagating along the background field lines. 
Thus, (apart from a weak non-linear wave-mode cou- 
pling) there is no preferential growth of current at the 
null point itself. Near the null, the current accumulates 
along the spine when the rotations disturb the fan plane, 
while rotations about the spine lead to currents in the 
location of the fan near the null. Such diffuse currents, 
extending along the separatrix surfaces, have been pre- 
dicted in incompressible analytic models 5, 6], although 
the form we find for the current sheets is somewhat dif- 
ferent. The plasma flows and field line reconnection as- 
sociated with these currents is of a rotational nature, as 
predicted by [Hj]. 

By contrast, when the spine or fan of the null is sub- 
jected to a shearing perturbation, there is a strong cur- 
rent growth which is localised at the null point itself. 
This is achieved by propagation of the disturbance across 
the background magnetic field lines, suggesting that the 
dominant wave mode associated with the disturbance is 
a fast mode. This behaviour agrees with the results of 
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FIG. 11: Evolution of two initially connected sets of field lines when the fan is sheared. One set are traced from plasma elements 
which remain always in the ideal region far up the spine (black, red online) and is advected across the spine. The other (grey, 
green online, traced from plasma elements initially within the black loop) are traced from plasma elements near the z — plane 
and are transported across the fan. Parameters as in Figure [5] but bo = 3, r\ constant = 0.002. The dashed field lines in the 
first image indicate the initial location of the spine and fan. Times for the images are t = 0, 0.8, 1.6 



(29| . and is consistent also with the work of [22| , who 
found singular current sheets to result at a 3D null point 
in an ideal line-tied relaxation when a shear was applied 
to spine or fan. The current which develops at the null 
is dominated by the component in the direction perpen- 
dicular to the plane of shear. The resulting plasma flow 
has a stagnation-point structure. It is in general very 
difficult to say from such simulations definitively that 
the plasma actually crosses the spine and fan. In the 
ideal region the flow advects field lines ideally, however 
it may still cross the spine or fan, if the spine and fan 
field lines change their identities (where now a field line 
is identified by plasma elements which lie on it in the 
ideal region) in time via magnetic reconnection (as in 
steady-state 2D reconnection) . Despite this difficulty, by 
comparison with kinematic models [25j , the non-zero in- 
tegrated parallel electric field which develops along fan 
field lines is a strong indicator of flux transfer across at 
least the fan. The magnetic field line evolution shows a 
reconnection of field lines through/around the spine and 
across the fan plane, as in the 'spine reconnection' and 
'fan reconnection' of [28] . The transport of magnetic flux 
across the fan plane is a particularly important property, 
as it is a separatrix surface of the magnetic field. The be- 
haviour of shearing perturbations is further in agreement 
with the kinematic predictions of [25| , 

Finally, we investigated the existence and behaviour 
of parallel electric fields within the simulations. In each 
case investigated, in order to halt the indefinite growth 
of J, resistive effects become important, as evidenced by 
the development of E\\ . In the case of rotational pertur- 



bations, the spatial profile of E\\ (as with J), is spread 
out either along the spine or the fan of the null. For 
shear perturbations, however, Eu tends to be strongly 
localised at the null point itself. Such parallel electric 
fields can be shown to give a physically meaningful mea- 
sure of the reconnection rate at a 3D null 0, l25l ] (though 
direct analogy with these kinematic results is not com- 
pletely straightforward here). Finally, the dependence of 
the behaviour of the system, and the scaling of the peak 
currents and reconnection rate, on resistivity, plasma-/? 
and the isotropy of the background null point are all im- 
portant issues, as is the energetics of the reconnection 
process. These will be left for exploration in a future 
paper. 
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